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Abstract 



^ . We consider the effects of uncertainties in nuclear reaction rates on the cosmological constraints on 

CO i the decays of unstable particles during or after Big-Bang nucleosynthesis (BBN). We identify the 

nuclear reactions due to non-thermal hadrons that are the most important in perturbing standard 
BBN, then quantify the uncertainties in these reactions and in the resulting light-element abun- 
f"~-« ■ dances. These results also indicate the key nuclear processes for which improved cross section data 

would allow different light-element abundances to be determined more accurately, thereby making 
possible more precise probes of BBN and evaluations of the cosmological constraints on unstable 
particles. Applying this analysis to models with unstable gravitinos decaying into neutralinos, we 
L^ calculate the likelihood function for the light-element abundances measured currently, taking into 

;_j ■ account the current experimental errors in the determinations of the relevant nuclear reaction rates. 

We find a region of the gravitino mass and abundance in which the abundances of deuterium, ^He 
and ^Li may be fit with x^ = 5.5, compared with x^ = 31.7 if the effects of gravitino decays are 
unimportant. The best-fit solution is improved to x^ ~ 2.0 when the lithium abundance is taken 
from globular cluster data. Some such re-evaluation of the observed light-element abundances 
and/or nuclear reaction rates would be needed if this region of gravitino parameters is to provide 
a complete solution to the cosmological ^Li problem. 



1 Introduction 

Late-decaying massive particles are generic features of plausible extensions of the Standard 
Model, such as supersymmetry. Cosmological constraints on such models are imposed by the 
observed astrophysical abundances of the light elements [1]- [23], which differ little from those 
calculated in standard Big-Bang Nucleosynthesis (BBN) [23H28] . An accurate calculation of 
the constraints imposed by astrophysical observations on the abundance of a late-decaying 
massive particle requires taking into account not only the uncertainties in the astrophysical 
observations but also the uncertainties in the nuclear reaction rates that contribute to the 
production of light elements in both standard and modified BBN scenarios. 

We report in this paper on a study of the effects on the astrophysical abundances of 
the light elements deuterium, ^He, ^He, ^Li and ^Li of the uncertainties in the rates for 36 
different nuclear reactions. Our central result can be expressed as a 36 x 5 matrix that may 
be used, e.g., to calculate the cumulative uncertainties in the BBN light-element abundances 
induced by the uncertainties in any given set of nuclear data, and to estimate the changes in 
the calculated abundances that would be induced by any updates of the measurements of the 
nuclear reaction rates. In particular, our analysis shows which uncertainties in reaction rates 
have the greatest impact on the light-element abundances calculated in standard BBN. Our 
analysis therefore bears upon the apparent ^Li problem for standard BBN [28j, and one of 
our main interests is in the application of our analysis to the constraints on the modifications 
to standard BBN that occur in models with late-decaying particles. 

Such particles appear, for example, in supersymmetric models with gravity-mediated 
supersymmetry breaking. If the gravitino is the lightest supersymmetric particle (LSP), 
then the next-to-lightest supersymmetric particle (NLSP) is relatively long-lived [29H3T] . 
Alternatively, if the lightest neutralino x is the LSP, then the gravitino is long-lived. In this 
paper we revisit the second possibility, as an illustration of the incorporation of the effects 
of the current uncertainties in the relevant nuclear reaction rates. 

Neglecting these uncertainties, we analyzed previously [32] the constrained minimal su- 
persymmetric extension of the Standard Model (CMSSM), in which the supersymmetry- 
breaking masses for spartners of Standard Model particles are assumed to be universal and 
the gravitino is assumed to be more massive [33]. We found that there were strips in rep- 
resentative gravitino mass vs. abundance (m3/2, C3/2) planes where the discrepancy between 
the measured ^Li abundance and that calculated in standard BBN could be reduced by the 
effects of late-decaying gravitinos without destroying the successful BBN predictions for the 
other light elements, particularly deuterium. 

In this paper, as an application of our general analysis of the implications of uncertainties 
in nuclear reaction rates, we include them in a re-evaluation of this possible supersymmetric 
solution to the ^Li problem. We re-calculate the global likelihood function x^ i^i the same 
representative (m3/2,C3/2) planes, now including the uncertainties in the measured abun- 
dances as well as the nuclear reaction rates. We confirm that x^ is indeed minimized along 
the strips found previously, being reduced typically by ~ 26 units. However, the quality 
of the best global fit to the light-element abundances is still not very good: x^ = 5.5 (for 
one effective degree of freedom - the three abundance measurements minus two fit param- 



eters). Thus, if late-decaying gravitinos are to solve the ^Li problem, they will need some 
help, either from changes in some reaction rates outside the uncertainties currently stated 
and/or changes in the measured deuterium and/or ^Li abundance, and we discuss some such 
possibilities. 

2 Principal Nuclear Reaction Rates 

We first discuss the nuclear reactions included in our analysis. The principle application of 
their non-thermal rates in the context of nucleosynthesis with late-decaying particles was 
discussed at length in [32], so here we review only some important details relevant to the 
current analysis. 

As described in [32j, non-thermal hadrons h G (p, n) are injected into the cosmological 
baryon/photon plasma by the decay of a heavy particle such as the NLSP, and then inter- 
act with the cosmic medium through which they travel. Non-thermal particle propagation 
is determined by competition among the various interactions that lead to particle losses- 
continuous energy losses as well as elastic and inelastic collisions. These loss processes are 
always rapid compared to the cosmic expansion rate. Thus, to a good approximation the 
non-thermal particle spectra are set by an equilibrium between injection and losses. These 
propagated, equilibrium spectra then determine the rates of non-thermal interactions with 
light nuclei via convolution with the relevant cross sections. That is, for the process hb ^ i 
of a non-thermal hadron interacting with a thermal background species b to produce light 
element i, the interaction rate per target b is Thh^i = J Nh{E) v ahb^i^E) dE, where Nh{E) 
is the spectrum of non-thermal h having kinetic energy E, and (Jhb^i is the cross section for 
the process at hand. The non-thermal processes considered here are listed in Table [H 

The uncertainty in the reaction rate Thi,^^, due to cross-section errors 6ahb^e, is 

SThb^, = j Nu{E) V 5aub^,{E) dE (1) 

= ^hb^l ^hb^i (2) 

where ehb^i = S^hb-^e/^hb^e characterizes the fractional error in the rate. The propagated 
non-thermal spectra Nh{E) generally increase to a peak at -E ~ few GeV, i.e., at energies far 
above thermal energies, the Gamow peak, and any reaction threshold. Non-thermal rates, 
unlike thermal rates, are sensitive to cross section behaviors over much larger ranges of 
energies. The cross sections typically grow rapidly above threshold, and then in some cases 
(fusion processes) drop strongly above E ~ few x 10 MeV, or in other cases (spallation 
processes) remain nearly constant or drop slowly at high energies. We should expect the 
uncertainties in non-thermal rates often to be larger than the typical uncertainties in the 
thermal rates, which are sensitive to a much narrower range of energies around the Gamow 
peak, typically ~ 0.1 - 0.3 MeV. 

In principle, the rates Thb^e have additional uncertainties due to those in the non-thermal 
spectra Nh, which would in general need to be added to the cross-section errors. Since the 
Nh are determined by sources and losses, they reflect the uncertainties in these processes. 



In practice, the dominant losses are typically electromagnetic energy losses to the plasma, 
the rates for which are relatively well known. Even in the regimes where the losses are 
dominated by scattering, the relevant cross-section errors are better known than typical 
reactions involving light elements. The source (e.g., NLSP decay) spectra are also well- 
determined by supersymmetry and Standard-Model physics. Consequently, the errors in 
Nh should be relatively small, and thus the uncertainties in the reaction rates should be 
dominated by the errors in the light-element cross sections that we have highlighted. 

We have estimated uncertainties for the non-thermal reactions by comparing nominal 
cross-section fits with experimental measurements. The fitting functions (j{E) typically 
provide good or excellent fits to the data. However, the data themselves are often sparse 
over the large energy ranges of interest. Unfortunately, this paucity of data is particularly 
acute for the spallation reactions /i^He — )■ (/i, ^A, ^A) + ■ ■ ■ , which are among the most 
important, as we shall see. In each case, we estimate conservatively the typical fractional 
size e = 6a /a of the experimental error bars over the energies where the cross section is 
substantial (i.e., near maxima for strongly-peaked cross sections, and out to ~ few GeV for 
flat cross sections). In the following section, we determine which of these reactions have the 
most important impacts on the light-element abundances, and we report the uncertainties 
for those in Table [H 

For a sufficiently large abundance of gravitinos, the standard BBN predictions are mod- 
ified, and the resulting light-element abundances need to be compared with observational 
determinations. In [32], we used the abundances (or abundance ratios) of D, ^He, ^Li, ^He/D, 
and ^Li/^Li to determine the allowed regions of parameter space defined by the gravitino 
mass, m3/2, the gaugino mass, mi/2, the ratio of Higgs vacuum expectation values, tan/3, 
and the gravitino abundance, ^3/2) characterized by 

7713/2^3/2 
C3/2 = = "^3/2^3/2??, (3) 

where 773/2 is the gravitino number density, Y3/2 = 773/2/^^, and rj = 6.19 x 10"^° is the 
baryon-to-photon ratio from WMAP year 7 [38]. For our present x^ analysis, we restrict our 
attention to the elements that have definite observational abundances with which we can 
make a comparison, namely the following. 

D/H: We use the deuterium abundance as determined in several high-redshift quasar 



absorption systems, which have a weighted mean abundance 

X 10"^ (4) 




where the uncertainty includes a scale factor of 1.7 due to the dispersion found in these 
observations. Since the D/H ratio shows considerable scatter, it is likely that systematic 
errors dominate the uncertainties. In this case it may be more appropriate to derive the 
uncertainty using sample variance (see e.g. [21]) which gives a more conservative range D/H 
= (2.82±0.53) X 10~^. We comment further on this below. The standard BBN result for D/H 
at the WMAP value for t] is (2.52 ± 0.17) x 10~^, showing potentially a slight discrepancy 
with the observed value, unless one adopts the larger uncertainty. 



Table 1: Nuclear reactions of non-thermal particles, including the most important of the 
estimated uncertainties in the cross sections. 



Code 


Reaction 


Uncertainty e 


Reference 


1 


p^He -^ d^Re 




Meyer [3^ 


2 


p'^B.e — )■ np^B.e 


20% 


Meyer 34j 


3 


p'^He — >■ ddp 


40% 


Meyer 34| 


4 


p'^Re — > dnpp 


40% 


Meyer 34 


5 


d^Re -^ *^Li7 




Mohr 3bJ 


6 


t^He ^ ^Lin 


20% 


Cyburt et al. Til 


7 


3He^He ^ ^Up 


20% 


Cyburt et al. l4j 


8 


t^He -> ^Li7 




Cyburt [27] 


9 


^He^He -^ ^Be7 




Cyburt and Davids j36 


10 


/Li -^ 3He^He 




Cyburt et al. 14 


11 


n^Li ^ t^He 




Cyburt et al. [l4j 


12 


pn — )■ (i7 




Ando, Cyburt, Hong, and Hyun 3Vj 


13 


pd — )■ ^He7 




Cyburt et al. [H] 


14 


pt -)■ ra^He 




Cyburt T\ 


15 


p^Li -> ^Be7 




Cyburt et al. |14| 


16 


p^Li -^ ®Be7 




Cyburt et al. 14| 


17 


p^Be ^ SB7 




Cyburt et al. 32] 


18 


np — )■ d'y 




Ando, Cyburt, Hong, and Hyun [37] 


19 


nd —7- t7 




Cyburt et al. 14 


20 


n^He -^ dt 




Meyer [3i] 


21 


n^He — >■ npt 


20% 


Meyer 34] 


22 


n^He — )■ ddn 


40% 


Meyer 34] 


23 


n^He — )■ dnnp 


40% 


Meyer 34] 


24 


nPLi -)■ ''Li7 




Cyburt et al. [13] 


25 


n (thermal) 






26 


n'^'Be -)■ p^Li 




Cyburt et al. |14J 


27 


n^Be ^ 4He4He 




Cyburt et al. 32] 


28 


p^Li ^ ^He^He 




Cyburt et al. 14] 


29 


mr~^ — )■ p7r° 




Meyer 34] 


30 


pvr" — )■ rnr^ 




Meyer 34] 


31 


p^He — )• ppt 


20% 


Meyer 34] 


32 


n'^He — )■ rm^He 


20% 


Meyer |34] 


33 


n^He — )■ nnnpp 




Meyer 34] 


34 


p^He — )■ nnppp 




Meyer 34] 


35 


p^He ^ A^^HeTT 




Meyer 34] 


36 


n^He -^ A^^HeTT 




Mever 341 



^He : The ^He abundance is determined from observations of extragalactic H II regions. 
These abundance determinations are known to suffer from large systematic uncertainties [16] . 



A recent analysis found [47] 

yp = 0.256 ±0.011, (5) 

and a similar central value was found in 08] . The standard BBN result for Yp at the WMAP 
value for rj is 0.2487 ± 0.0002, which is consistent with observations, given the error in ([5]). 
^Li/H: The ^Li abundance is derived from observations of low-metallicity halo dwarf 
stars. Some > 100 such stars show a plateau 09] in (elemental) lithium versus metallicity, 
with a small scatter consistent with observational uncertainties. An analysis [SD] of field halo 
stars gives a plateau abundance of 

w) = (l-23l^i6) X 10-1°, (6) 

^ / halo* 

where the errors include both statistical and systematic uncertainties. As in the case of ^He, 
the errors are dominated by systematic uncertainties. For example, the lithium abundance 
in several globular clusters, tends to be somewhat higher [^TH5B] . and we make some com- 
parisons below to the result found in [56] of ^Li/H = (2.34 ± 0.05) x lO-^". However, the 
standard BBN result for ^Li/H at the WMAP value for rj is (5.12+|^:^2) x 10"^°, which differs 
significantly from the observed value, hence the ^Li problem [28]. Note that the central 
values for the BBN abundances used here differ slightly from those in fIE\ , primarily due to 
the small shift in rj as reported in [58] . 

Recently, there have been several analyses which indicate that the ^Li abundance at 
low metallicity falls below the typical plateau value and/or shows a significant amount of 
dispersion [57I - I6T] . While these observations apparently provide the first indications of Li 
depletion in metal-poor stars, it would appear that it is operative only at extremely low 
metallicity, [Fe/H] < —3, whatever the particular depletion mechanism may be, whether in 
the star or in the medium prior to the star's formation. There is no observational evidence 
of any depletion at higher metallicity (—3 < [Fe/H] < —1.5) from the standard BBN result 
to the plateau value ^59ii60J, in contrast to the claim of [61j [J. 

To obtain our x^ distribution, we combine the standard BBN uncertainties with the 
observational errors in quadrature. In the case of ^Li, where the reported errors are uneven, 
we use the upper error bar on the observation, and the lower error bar on the theory, as we 
are interested in the region between these two central values. Correspondingly, the likelihood 
function that we calculate is 

2 ^ f Yp~0.25Q y / g-2.82x lO'^ V / ^ - 1.23 x lO'^" . ,-^ ^ 
^ ^ V 0.011 / ^ I 0.27 X 10-5 / ^ 0.71 X 10-10 I + Z^ ^^ ' ^'> 



iJt was argued in "pT that there are two plateau values corresponding to [Fe/H] above and below -2.5. 
However, the evidence for this assertion is not convincing, as these data can be fit equally well with a linear 
increase in logLi vs. [Fe/H] as in (SDHSHUHOIISS] ■ This would point to a lower primordial Li abundance and a 
more severe problem with respect to standard BBN predictions. 



where the Sj are the contributions to the total x^ due to the nuisance parameters associated 
with varying one or more of the rates hsted in Table [H Standard BBN has a large total 
X^ = 31.7, primarily due to the discrepancy in ^Li. There is a contribution of Ax^ ~ 30 
from the ^Li abundance, Ax^ ~ 1.2 from the D/H abundance, and a smaller contribution 
from ''He, corresponding to a ~ 5 — cr discrepancy overall q 

Our treatments of the hadronic and electromagnetic components of the showers induced 
by heavy-particle decays follow those in [32] . Also, we follow the calculations of decay 
branching ratios and particle spectra described in |32]. The only differences here are in the 
nuclear reaction rates and their uncertainties that were discussed above. 

We display in Fig. [1] the effects on the abundances of the light elements deuterium, 
^He, ^He, ^Li and ^Li of the decays of a generic metastable particle X with lifetime tx G 
(1, 10^°) sec. For illustration, we assume the decay spectra calculated in [32] for the choice 
(mi/2,m3/2, tan/3) = (300 GeV, 500 GeV, 10), in which case the proton branching ratio 
Bp ca 0.2 and the electromagnetic branching rate is -Bem"^3/2 = 115 GeV. In this figure we 
assume the nominal central values of the nuclear reaction rates discussed in the text, and 
this figure may be compared directly with Fig. 6 of [32j. The main differences are in the 
upper left panel, where the region where the deuterium abundance lies within the favoured 
range is now pushed to values of (x that are lower by a factor of about 2 when tx < 10® sec 
as compared with the results of [32], and in the lower middle panel, where the region of 
acceptable ^Li abundance extends to lower (x when tx ~ 10^ sec. Both these effects are 
due to the inclusion of the reactions n^He — )■ nn^He and p^He — ?• ppt, and have the effect of 
pushing the location of a possible 'solution' of the ^Li problem also to lower (x- 

We show in Fig. [21 one generic (m3/2,C3/2) plane, also without the inclusion of uncer- 
tainties in the non-thermal rates in Table [H This plot is based on a specific CMSSM point 
(benchmark C of [S3]) with mi/2 = 400 GeV, Aq = 0, and tan/3 = 10. The universal scalar 
mass is set to niQ = 90 GeV to get the correct WMAP density for dark matter. The lightest 
neutralino mass is about 165 GeV for this point, and for gravitino masses larger than this 
we have neutralino dark matter with an unstable massive gravitino. As in Fig. [T|, the region 
where the deuterium abundance lies within the favoured range is now also pushed to lower 
^3/2 when m3/2 ^ 2 TeV as compared to the results in |32j, and the region of acceptable 
^Li abundance extends to lower (^3/2 when 7713/2 is between 2-3 TeV. In the lower right 
panel, we see marginal compatibility between the ^Li constraint (light blue) and the other 
constraints for 7713/2 ^ 3 TeV. This region will be the focus of our discussion in the following 
X^ analysis. 

3 Incorporation of Uncertainties 

Of the 36 interactions we study, there are just 10 whose uncertainties induce non-negligible 
uncertainties in the light-element abundances, namely the reactions 2, 3, 4, 6, 7, 21, 22, 23, 
31, and 32 in Table [TJ Their uncertainties are not important for the ^He abundance Yp, but 



^We find that x^ ~ 21.8 even when the globular cluster value of ^Li/H is used, corresponding to a 4 — ct 
effect. 
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Figure 1: Plots of the effects on the light- element abundances of the decays of a generic 
metastahle particle X with a lifetime tx G (1, 10"'^°) sec, using the decay spectra calculated 
for (mi/2,m3/2, tan/3) = (300 GeV, 500 GeV, 10), in which case Bp ^ 0.2 and the electro- 
magnetic branching rate is -Bem^™3/2 = 115 GeV. The X abundance before decay is given 
by Cx = f^xiT-x/iT-y- The white regions in each panel are those allowed at face value by 
the ranges of the light-element abundances reviewed in Section\^ whilst the yellow, red and 
magenta regions correspond to progressively larger deviations from the central values of the 
abundances. 



are potentially important for the deuterium, '^He, ^Li and ^Li abundances. 

To explore the effect of reaction uncertainties, we quantify the reaction sensitivity as 
follows. Using the set of unperturbed non-thermal reaction rates {r°}, we find the unper- 
turbed abundances of light elements: examples of these results appear in Figures [1] and [2l 
For a given light element i, we call the unperturbed abundance y^. Then, for any single 
reaction j, we consider changes in the rate by a factor 1 + e: F' = (1 + e)F°, leaving all other 
non-thermal (and thermal) rates unchanged. We evaluate the new resulting light-element 
abundances for a wide range of values for e including both positive and negative values and 
write the new, perturbed i abundance as y'llrxnj = Vi + Syilrxnj- In this way, we were able to 
identify the 10 reactions listed above as potentially playing an important role in altering the 
light-element abundances. Our final results are based on a Gaussian distribution of rates 
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Figure 2: The effects of the decays of a gravitino with variable mass 771^/2 on the different 
light-element abundances for a specific point (benchmark C of l6^ ) with mi/2 = 400 GeV 
on the WMAP coannihilation strip for a CMSSM scenario with tan/3 = 10, Aq = 0. As in 
the previous figure, the white regions in each panel are those allowed at face value by the 
light-element abundances reviewed in Section\^ and the yellow, red, and magenta regions 
correspond to progressively larger deviations from the central values of the abundances. 



with widths give by the values of e chosen according to the uncertainty estimates in Table [U 
We display in Figs. |3l HI O and [7] the effects of the uncertainties in these reaction rates 
on the abundances of each of the key elements among deuterium, ^He, ^Li and ^Li, each in 
a (^3/2; C3/2) plane. We concentrate here on benchmark point C, as the effect of perturbing 
the interactions is qualitatively similar for the other benchmark points we consider below. 

For example, in Fig. E] we show the effect of the p'^He — )■ rap^He (reaction 2) in Table [T] 
on the abundances of D/H (left) and ^He/H (right). In this case the effect on ^Li and ^Li is 
negligible. We plot contours showing 
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the relative change in the light-element abundance when rate j is perturbed by a factor 
(1 + e). For reaction 2, we estimate a 20% uncertainty in the rate and, as one can see, the 
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Figure 3: The effects in the (^713/2, C3/2) plane of the 20% uncertainty in the rate for the 
reaction 2 (p^He — ?■ np^Yie) on the abundances of deuterium (left) and^Yie (right). Contours 
show the relative changes in the light-element abundances. 



effect on D/H is always less than 4% and occurs at very high gravitino abundances. The 
effect on ^He is also relatively small, but extends over a larger portion of the parameter 
space. The effect of reaction 31( p^He — )■ ppt) is qualitatively similar to the one shown in 
Fig. [31 

For reactions 3 and 4, corresponding to p^He — )■ ddp and p'^He — > dnpp respectively, we es- 
timate an uncertainty of 40%. However, even with the larger uncertainty, the relative change 
in D/H for both rates is still less than 4%, extends down to lower C3/2 ~ 10~^° GeV, and the 
^He abundance variation is even smaller. Accordingly, we do not show these examples. 

In Fig. m the effects of reaction 6 corresponding to f^He — )■ ^Lira are displayed (the effects 
of reaction 7 corresponding to ^He^He — )■ ^Lip are similar but weaker by a factor of 2): we 
estimate 20% uncertainties for these reactions. We see that in this case, while there is some 
effect on the abundance of ^Li, the dominant effect of varying this rate is on ^Li, where 
changes can be as large as 12% for almost all the values of C3/2 shown when 7713/2 ~ 1 TeV. 
We find similar results for points E and L, whilst for point M (see below) similarly large 
changes in ^Li are centered around 7713/2 ~ 2 TeV. The effects on deuterium and ^He are 
negligible for reactions 6 and 7. 

The effects of reaction 21 (77^ He — )■ npt)., for which we also estimate an uncertainty of 
20%, on all four light elements are shown in Fig. [51 The possible effect on ^Li is largest, 
amounting possibly to a reduction in the ^Li abundance by up to 6% in a diagonal region 
extending from 7^3/2 > 3 TeV. (For benchmark point M, the reduction in the ^Li abundance 
occurs at 777,3/2 > 4 TeV.) Reaction 32 (77"^ He — )■ r777^He) shows effects somewhat similar 
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Figure 4: Similar to Fig. 0, /or i^/ie reaction 6 ft^He — )■ ^Linj, m t/izs case showing the effect 
on '^U/H (left) and ^U/H (right). 

(though in general a bit weaker) to those seen in Fig. [5], and is not shown separately. 

We note that, although the relative shift in the abundance is small for a 20% variation 
in the rate, the ^Li abundance in this region (the diagonal strip where the abundance is 
decreased by 4-6% in Fig. [5]) is already significantly reduced when using the (unperturbed) 
non-thermal rates. In this region, several rates (principally reactions 21, 23, and 32) com- 
bine to lower the ^Li abundance when the gravitino abundance is sufficiently large. However, 
subsequent variations in the rates do not make any further significant changes in the abun- 
dances. To help better understand this point quantitatively, we show in Fig. El the ^Li 
abundance as a function of e for rate 21 (other rates have ei^2i = 0) in the upper panel 
and as a function of €21 = 623 = £32 in the lower panel. As one can see, particularly in 
the latter case, when €21,23,32 = —1 and these rates are shut off entirely, the abundance of 
^Li is 4.2 xlO~^° (the 20% decrease in ^Li is due to the remaining non-thermal reactions). 
Furthermore, coherent variations in these rates of 20-40% make relatively small changes in 
the abundances as refiected in Fig. [5], and the effects of random variations in the rates would 
clearly be smaller still. 

Whilst reaction 23 (n^He — )■ dnnp) shows smaller variations in ^Li, the uncertainty (which 
we estimate at 40%) is larger. The effect on D/H is also pronounced, as seen in Fig. [71 The 
effect of reaction 22 (n^He — )■ ddn) is qualitatively similar but weaker. 

We conclude this Section by summarizing the main results of our propagation of non- 
thermal reaction rate errors into uncertainties in light-element abundances. We find that 
the *^He abundance is essentially unaffected by reaction rate errors. For ^Li, we find that no 
one reaction dominates the non-thermal perturbations, which in turn means that errors in 
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four light elements deuterium (upper left), '^He (upper right), '^Li (lower left) and ^Li (lower 
right). 
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any given rate only have a rather small (typically < 10%) effect on the ^Li/H abundance. 
Non-thermal deuterium production is also not entirely controlled by a single reaction, though 
the n^He — )■ dnnp reaction clearly stands out as the most important, and the resulting errors 
in D/H can go as high as > 20%. 

These results have important implications for our x^ analysis, which, as we will see, is 
dominated by ^Li/H and D/H. Since the ^Li/H nuclear uncertainties are small compared to 
the observational errors in the ^Li/H abundance, the latter dominates the lithium contribu- 
tion to x^- Conversely, the D/H non-thermal rate errors are significant in comparison to the 
observational errors, and thus will have an important effect on the x^ ^ind ultimately on the 
lithium problem. From this we infer that the reactions which most critically need improved 
nuclear data are those which are important for deuterium production. 

4 x^ Analyses of Benchmark CMSSM Scenarios 

To proceed with the x^ analysis, we use Eq. I^^ to calculate x^ for each point sampled in the 
(^^3/2,C3/2) plane. The reaction rates are treated as nuisance parameters and therefore, for 
each evaluation of x^? each non-thermal rate is chosen from a Gaussian distribution about 
the mean rate with the uncertainty specified in the previous Section. At each point and for 
each reaction considered, the difference between the rate chosen and its mean value, relative 
to the quoted uncertainty in the rate, determines the corresponding Si in Eq. (jTj). 

From the results of the analysis in the previous Section, it is clear that it will be sufficient 
to focus on the effects of reactions 21, 23, and 32. In principle, one could include all reactions 
in the x^ analysis as nuisance parameters. However, the inclusion of many more reactions 
would have only a marginal effect on lowering the x^ contribution from the abundances, 
while at the same time increasing x^ through sj. Since each rate typically increases x^ by 
roughly one unit, one would need to gain at least one unit from the effect of the uncertainty 
in the rate on the element abundances. Including the uncertainties of reactions beyond 21, 
23, and 32 with finite sample sizes will typically lead to a larger value of x^- 

In [22], we discussed the application of the BBN constraints to four benchmark CMSSM 
scenarios with specific values of the soft supersymmetry-breaking parameters mi/2,mQ and 
Aq, tan/?, and the Higgs mixing parameter /i, labelled C, E, L and M |63j. In each case, 
Afj = and /i > was chosen. The parameters corresponding to point C were given earlier. 
For points E, L, and M they are (mi/2,mo, tan/3) = (300 GeV,1615 GeV, 10), (460 GeV, 310 
GeV, 50), and (1840 GeV, 1400 GeV, 50) respectively. Variants of these CMSSM scenarios 
with a massive gravitino are characterized by the gravitino mass m3/2 and its abundance 
(3/2- The (m3/2,C3/2) planes for benchmark scenarios C, E, L and M are shown in Fig. [HI 
displaying x^ contours for the light-element abundances calculated incorporating the nuclear 
reaction rate uncertainties discussed above. 

In the limit of large ms/2 and/or small C3/2, the value of the x^ function approaches 
~ 31.7, the same value as in standard BBN. This large value of x^ is due primarily to the ^Li 
problem. We see that in each of the CMSSM scenarios in Fig. [H] there is a 'trough' of much 
lower x^ with a minimum at ~ 5.5, shown in each panel by a cross. We display contours 
of x^ = 6 and 9.2, corresponding to the 95 and 99% CLs for fitting to two parameters. 

13 




-9 



(N 



m 




J^» 


-10 


bU 




O 




U 


-11 




-12 




-13 



, ,, ^ , , , , 


1 o.oe/ 


,, , , 


\__.^^^ ' ' - 


\\-0.04 


1 J 




0.02 




|\ 


0.04^-- 


-^ j^ai : 


- 


-0.02 


k ^ 


- 


- 






^He/H - 





, 









bX) ' 

O 



10 



11 



-12 



-13 



12 3 4 

m3/2 (TeV) 



12 3 4 5 

m3/2 (TeV) 




"Li/H 



2 
^3/2 



ma/o (TeV) 



Figure 7: Similar to Fig. 0, /or it/ie reaction 23 (n'^He ^ dnnp). 



14 




2000 



3000 4000 

m3/2 (GeV) 



5000 2000 



3000 



4000 



■^3/2 (GeV) 

















-8- 


Point L 




































,•• * 


























































-9- 






^ 


<^ 


^ 




r. -10- 




,.-•;■''' -- 


/ 






- 
















J^J' 




„„.;•'' ''./^ 


^ 








60 


50 ....■;.- 


y""^ ^ 


^ 








o 














J -11- 


.;...-3. ^ 


_^--^^ 










-12- 












.' 
















-13- 















O 



2000 3000 4000 

■^3/2 (GeV) 



5000 



3000 



4000 
m3/2 (GeV) 



5000 




5000 
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the uncertainties in the nuclear rates discussed in the text. 
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Table 2: Results for the best-fit points for CMSSM benchmarks C, E, L and M. The second 
set of results for C and M correspond to the globular cluster value for primordial ^Li/H. The 
third and fourth entries for point C correspond to the higher adopted uncertainty for D/H 
in field stars and to the globular cluster ^Li abundances, respectively. 



7713/2 GeV Logio(C3/2/[GeV ) 


Yv 


D/H (x 10-5) ^Li/H (x 10-10) Y.si 


2 

X 


BBN 


0.2487 


2.52 5.12 


31.7 



C 4380 -9.69 0.2487 3.15 2.53 0.26 5.5 

E 4850 -9.27 0.2487 3.20 2.42 0.29 5.5 

L 4380 -9.69 0.2487 3.21 2.37 0.26 5.4 

M 4860 -10.29 0.2487 3.23 2.51 1.06 7.0 



c 

M 


4680 
4850 


-9.39 

-10.47 


0.2487 
0.2487 


3.06 
3.11 


2.85 
2.97 


0.08 
0.09 


2.0 

2.7 


C 


3900 


-10.05 


0.2487 


3.56 


1.81 


0.02 


2.8 


c 


4660 


-9.27 


0.2487 


3.20 


2.45 


0.16 


1.1 



Also shown are the higher x^ contours of 32 (corresponding to the BBN value) and 50. We 
see that the (^713/2, C3/2) planes are very similar for benchmarks C, E and L. The plane for 
benchmark M is somewhat different, and the minimum value of x^ is slightly higher. In 
Table [21 we show the various abundances and x^ contributions for each of the three light 
elements for the standard BBN result and our best-fit point for each of the four benchmark 
points. 

It is interesting to note the tension between D and ^Li. At each of the best fit points, there 
is a considerable reduction in ^Li, approaching the observational value. The minimum value 
X^ ~ 5.5 certainly amounts to a 'mitigation' of the ^Li problem, but not a 'solution', in the 
sense that since we are fitting two parameters (m3/2 and C3/2) and using 3 measurements, we 
have effectively only one degree of freedom and x^/d.o.f . ~ 5.5. However, this improvement 
in ^Li comes at the expense of D/H, which at this point begins to make a more significant 
contribution to the total x^- On the other hand, the "^He abundance Yp does not contribute 
significantly to the likelihood at any point in the parameter space. At the minimum, the 
deuterium abundance contributes Ax^ ~ 1.5, whereas the '^Li abundance contributes Ax^ ~ 
3.4. Thus the previous 4- or 5-cr ^Li problem is reduced to a < 2-cr problem. If this mitigation 
is to lead to a complete solution, one or more of the nuclear reaction rates and/or measured 
light-element abundances should lie outside its quoted uncertainty. 

As an example, in Fig. [Hlwe show the analogous results for the x^ likelihood, assuming 
the globular cluster value for ^Li/H. Results for this case for benchmark points C and M are 
also summarized in Table |2l We now see the appearance of contours for x^ = 4.6 and 2.3 
corresponding to 90 and 68 % CLs respectively. The best-fit x^ values drop considerably in 
this case, with values of 2.0 and 2.7 for points C and M respectively. Thus a massive (> 4 
TeV) gravitino can provide a potential solution of the lithium problem if globular cluster 
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Figure 9: As in Fig.l^ contours of the x^ function in the {1713/2, C3/2) planes for the benchmark 
CMSSM scenario C (left) and M (right), assuming the globular cluster value of ^Li/H. 



data is assumed to represent the primordial ^Li abundance. 

As discussed earlier, one may also consider the effect of increasing the size of the uncer- 
tainty in the mean D/H abundance. Using an observed abundance of (2.82 ± 0.53) x 10~^, 
we obtain the x^ contours seen in the left panel of Fig. [101 corresponding to point C In this 
case, we can obtain solutions with x^ = 2.8 and a best-fit point with a ^Li/H abundance of 
1.81 xlO~^° coming at the expense of a higher D/H abundance of 3.56 x 10~^. When the 
globular cluster value of ^Li/H is used together with the higher D/H uncertainty, we can 



even find a best-fit solution with x^ 
seen in the right panel of Fig. [101 



1.1: D/H = 3.20 x 10"^ and ^Li/H = 2.45 x 10" 
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as 



5 Summary and Conclusions 

We have presented in this paper an analysis of the modifications of the cosmological light- 
element abundances that would be induced by the late decays of massive particles, incorpo- 
rating for the first time the uncertainties in relevant nuclear reaction rates. We have analyzed 
the possible effects of the 36 different nuclear reactions shown in Table [H and identified three 
as the most important, namely n'^He — )■ npt, n^He — )■ dnnp and n^He — )■ nn'^He. 

It is well known that there is a problem with the cosmological abundance of ^Li in conven- 
tional BBN with no late-decaying particles, and a natural question is whether this problem 
could be mitigated by some suitable late-decaying particle. As an example of the possible 
applications of our uncertainty analysis, we have considered in this paper the late decays of 
massive gravitinos in various benchmark supersymmetric scenarios. It had been observed 
previously that there were regions of the parameter spaces of these scenarios, corresponding 
to ranges of ms/2 and C3/2, where the cosmological ^Li problem might indeed be mitigated. 
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mark CMSSM scenario C assuming (left) a greater uncertainty in the observed D/H abun- 
dance and (right) also assuming the globular cluster value of ^Li/H. 



and we have made a likelihood analysis of this possibility incorporating uncertainties in the 
nuclear reaction rates. 

We confirm that there are indeed regions of the (m3/2, C3/2) parameter planes in these 
scenarios where the global x^ function is reduced from its value ~ 31.7 in conventional 
BBN (~ 21.8 if the globular-cluster value for the ^Li abundance is adopted) to x^ ~ 5.5. 
This provides a very significant alleviation of the ^Li problem, reducing it from a 4- or 5-a 
problem to a < 2-cr issue. The fact that our best-fit points lie within the x^ = 6 contours in 
Fig. [8] implies that they have a goodness-of-fit slightly exceeding 5%, which is marginal for 
considering massive gravitino decay as a 'solution' to the cosmological *^Li problem. 

The fact that this prospective solution exists for several choices of supersymmetric sce- 
narios, with parameters that are relatively stable, suggests that it is a general feature of 
supersymmetric models. For this potential solution to be confirmed, one or more of the 
following should happen. 

1. There might be some refinement in measurements of the cosmological ^Li abundance 
leading to a shift in the central value and/or a change in the assigned uncertainty. 
As we have shown in Fig. [HI for example, if the globular-cluster estimate of the ^Li 
abundance is adopted (which would correspond to x^ ~ 21.8 in standard BBN), the 
decays of massive gravitinos could reduce x^ to ~ 2.0. 

2. Alternatively, it is possible that the rates for one or more nuclear reactions might lie 
outside the ranges favoured by the current measurements and their assigned uncertain- 
ties. As we have pointed out, the measurements of some of these non-thermal rates are 
sparse over the energy ranges of interest, and improved coverage is certainly possible 



and desirable. The highest-priority reactions for new cross-section measurements are 
n^He — !■ npt, n'^B.e — )• dnnp and n'^He — )■ nn^He, which we have shown to be the most 
relevant for this analysis. 

3. Finally, we should mention the possibility of some unidentified error in our analysis: 
we have given reasons why we think its uncertainties are smaller than those mentioned 
earlier in this paragraph, but we could be wrong. 

The supersymmetric possibility of 'solving' the cosmological ^Li problem is currently at a 
very intriguing stage. The decays of massive gravitinos are one possibility, but it is a generic 
feature of supersymmetric theories with gravity-mediated supersymmetry breaking that there 
is some late-decaying massive particle, and other possible candidates exist. We plan to 
return to some possibilities in a forthcoming paper. Clearly, this approach to 'solving' the 
cosmological ^Li problem would be given an enormous boost if experimental evidence were 
to emerge for supersymmetry, either at the LHC or in (in) direct searches for astrophysical 
dark matter. In this connection, we note that there are good prospects for discovering 
supersymmetry at the LHC in many benchmark scenarios, including points C, E and L 
discussed here [63], and that there are also promising prospects for (in) direct detection of 
supersymmetric dark matter in scenarios C, E and L [61]. If supersymmetry were to be 
discovered, the search for evidence of a possible metastable supersymmetric particle would 
assume high priority, and it would be an exciting challenge to correlate its possible roles in 
cosmology and in the laboratory. Then one might indeed be justified in claiming that the 
cosmological ^Li problem was solved. 
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